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ABSTRACT 



A model for GRMHD disk outflows with neutrino-driven mass ejection is developed, and 
employed to calculate the structure of the outflow in the sub-slow magnetosonic region and the 
mass loading of the outflow, under conditions anticipated in the central engines of gamma-ray 
bursts. The dependence of the mass flux on the conditions in the disk, on magnetic field 
geometry, and on other factors is carefully examined for a range of neutrino luminosities expected 
in hyperaccreting black holes. It is found that a fraction of up to a few percent of the neutrino 
luminosity produced inside the disk is deposited in the upper disk layers, resulting in a steep rise 
of the specific entropy and the consequent ejection of baryons along magnetic field lines. For the 
range of conditions explored the final value of the dimensionless entropy per baryon in the wind 
is typically below 100. The fraction of neutrino luminosity that is ultimately being converted 
to kinetic energy flux is shown to be a sensitive function of the effective neutrino temperature 
at the flow injection point, and the shape of magnetic field lines in the sub-slow region, but is 
practically independent of the strength of poloidal and toroidal magnetic fields. We conclude 
that magnetic launching of ultra-relativistic polar outflows from the innermost parts of the disk 
is in principle possible provided the neutrino luminosity is sufficiently low, L v ,$ 10 52 erg s^ 1 
or so. The conditions found to be optimal for the launching of an ultra-relativistic jet are also 
the conditions favorable for large neutron-to-proton ratio in the disk, suggesting that a large 
neutron excess in GRB jets, as often conjectured, may be possible. However, the outflow time 
appears to be comparable to the neutronization timescale, implying that the electron fraction 
should evolve during the initial acceleration phase. Further analysis is required to determine the 
composition profile in the wind. 



Subject headings: accretion, accretion disks - gamma rays: bursts - MHD - relativity - nuclear 
reactions, nucleosynthesis, abundances 



1. Introduction 

The central engines that power gamma-ray bursts (GRBs) are thought to consist of a newly formed 
stellar mass black hole surrounded by a hot, dense, magnetized disk. The prompt and afterglow emissions 
are inferred to be produced in ultra-relativistic outflows that are launched by the central engine and 
accelerate to Lorentz factors 7^ > 100. The free energy source of those relativistic, gamma-ray emitting 
jets can be either the gravitational potential energy of accreted matter or the spin energy of a Kerr black 
hole. Whether due to dissipation of some fraction of the spin energy of the hole (van Putten & Levinson 
2003; Levinson 2005), or owing to hyper-accretion rates (Popham et al. 1999; Pruet et al. 2003), the disk 
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surrounding the black hole is expected to be sustained at MeV temperatures, and its midplane density 
is expected to be high, in excess of 10 10 gr cm~ 3 . Under such conditions the disk is optically thick to 
electromagnetic radiation and the dominant energy loss mechanism is neutrino emission. Additionally, the 
nuclear composition in the inner parts of the disk is likely to be neutron rich (Bcloborodov, 2003; Pruct 
et al. 2003); the neutron-to-proton ratio can in principle approach 30 if the disk material is cold enough 
(T £ 2 MeV) and dense enough (p £ 10 11 gr cm~ 3 at the disk midplane; see e.g., fig. 1 in Bcloborodov 
[2003]). The prodigious neutrino luminosity is expected to drive a powerful wind from the disk. This wind 
is most likely baryon rich and expands at sub or mildly relativistic speeds except, perhaps, inside a core 
containing the putative baryon poor GRB jet. The surrounding baryon rich wind may play an important 
role in the collimation of the central jet (Levinson & Eichler, 2000; Rosswog & Ramirez- Ruiz, 2003), and 
its presence may also lead to some observational consequences (e.g., a supernova-like event). The nuclear 
composition of the disk outflow, particularly the central GRB jet, is an issue of considerable interest. If, 
as often assumed, the GRB jet picks up nuclei from the disk (e.g., Derishev, Kocharovsky & Kocharovsky 
1999; Beloborodov 2003) then it may contain matter with a large neutron excess that is likely to affect 
its dynamics (Fuller et al. 2000; Vlahakis, Peng & Konigl 2003; Rossi et al. 2005), as well as some of 
the characteristics of the afterglow emission (e.g., Derishev, Kocharovsky & Kocharovsky 1999; Bachall & 
Mcszaros 2000). In addition, the disk wind may be a potential site for nucleosynthesis (Pruet et al. 2004, 
2005). The nuclear composition of the disk outflow should depend on the conditions at the disk surface and 
the details of mass ejection process, particularly on the ratio of outflow time and neutronization timescale. 

The gross features of the physical picture outlined above are supported by recent numerical simulations 
(e.g., Proga et al. 2003; McKinney, 2005a; and references therein), however, the detailed structure and 
properties of the polar outflow are only marginally resolved. In particular, the origin of the fireball - 
whether emanates from the inner regions of the disk or produced via a BZ mechanism, the nature of the 
baryon loading process, and the nuclear composition of the outflow are yet open issues. Moreover, as will 
be shown below the mass ejection process is rather sensitive to the conditions at the flow injection point 
and the physical processes involved, which makes it difficult to study using global numerical simulations. 

In order that a wind be accelerated hydrodynamically to high Lorentz factors, the entropy per 
baryon at the wind injection region must be extremely high. Purely hydrodynamic fireballs demand 
s/ks ^ 10 5 (joo/'200){kTd/2MeV)~ 1 where is the terminal Lorentz factor of the fireball and Td is 
the temperature at the base of the outflow. Although the specific entropy of the fireball should rise as 
it accelerates, owing to absorption of neutrinos escaping from the dense layers of the disk, such extreme 
values are unlikely to be achieved. The requirements on the baryon loading can be alleviated if the outflow 
is magnetically dominated. In that case, the disk must support an ordered magnetic field with strength 
in excess of a few times 10 14 G to account for the observed GRB energies (e.g., Levinson & Eichler 1993; 
Meszaros & Rees, 1997). Such strong fields conceivably exist in the dense disks of hyperaccreting black holes 
(e.g., McKinney, 2005b). Since the energy is extracted magnetically, acceleration of the outflow to high 
Lorentz factors does not require, in principle, high specific entropy, but rather a high magnetic energy per 
baryon Sb- In the ideal MHD case Sb > 7oo ~ 200 and if the fireball accelerates along magnetic field lines 
that are anchored to the disk, the question then arises as to how baryon overloading of magnetic field lines 
can be avoided. Vlahakis, Peng & Konigl (2003) proposed that such high Eb may not in fact be necessary if 
the ratio of baryon-to-magnctic flux is not conserved on magnetic flux surfaces. Specifically, they considered 
a situation in which a neutron rich outflow, expelled along magnetic field lines that emanate from accreted 
disk material in the vicinity of the innermost stable circular orbit, is accelerated magnetically to a Lorentz 
factor of order 10, at which point the neutrons decouple. The remaining protons then continue to accelerate 
magnetically to <~ 200. With n/p ratio of about 30, £b ~ 10 or so is sufficient before decoupling (but 
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not smaller, otherwise the neutrons wont decouple at all). However, such a high n/p ratio requires optimal 
conditions in the disk, and essentially no evolution of the electron fraction Y e during the acceleration phase 
prior to decoupling, which is questionable. 

The issue of mass loading has not been addressed in most of the previous work on magnetized, 
relativistic disk outflows. The common approach, following the pioneering work of Blandford & Payne 
(1982), is to seek self-similar solutions of the trans-field equation (Li et al., 1992; Contopoulos 1994; Vlahakis 
& Koniglc 2003). While such a treatment considerably simplifies the analysis, it precludes the incorporation 
of gravity (in the relativistic case) as well as neutrino heating in the flow equations, and in addition does 
not allow matching of the self-similar outflow solution to a Keplerian disk (Li et al. 1992). The immediate 
implication is that the self-similar outflow solutions cannot be extended down to the sub-slow magnetosonic 
region. On the other hand, the mass-to-magnetic flux ratio is determined by the regularity condition at 
the slow magnetosonic point, namely by the requirement that the outflow must pass smoothly through 
this point. This means that solutions of the flow equations in the sub-slow magnetosonic region must be 
obtained in order for the mass flux to be calculated. Since the self-similar treatment is inapplicable in this 
region, one must seek a different approach. In this paper we construct a model for GRMHD disk outflows 
with neutrino-driven mass ejection, and employ it to calculate the structure of the flow in the sub-slow 
region. Using the regularity condition at the slow magnetosonic point we carefully examine how the mass 
flux depends on the conditions at the injection point, on geometry of magnetic surfaces, and on other wind 
quantities. In §2 we present the basic equations for a GRMHD wind with external energy and momentum 
sources. In §3 we derive analytic expressions for the flow quantities in terms of the location and temperature 
at the slow point by employing the regularity condition, and identify some systematic trends. In §4 we 
present numerical solutions. Our approach is to integrate the flow equations assuming the magnetic field 
geometry is given. Since the slow point is located close to the disk midplane, our choice of a particular field 
geometry reflects in part a choice of boundary conditions. We explore a range of field geometries in order 
to elucidate the dependence of the mass flux on the shape of the field lines. 

Neutrino-driven winds have been discussed previously in other contexts (e.g., Duncan et al. 1986; 
Levinson & Eichler 1993; Woosley et al. 1994; Witti et al. 1994; Qian & Woosley 1996; Thompson et al. 
2001). However, most previous models assume purely hydrodynamic, spherically symmetric winds, which 
are less relevant to the problem at hand. A simplified model of hydrodynamic disk outflow is presented in 
Pruet et al. (2004), mainly in connection with nucleosynthesis. In this work we generalize those previous 
treatments to incorporate GRMHD effects and disk geometries. 

2. GRMHD Neutrino- Assisted Wind Model 

We consider a MHD wind expelled from the surface of a hot, magnetized disk surrounding a 
non-rotating black hole. The range of conditions in the disk is envisioned to be similar to that computed 
by Popham et. al. (1999) for hyperaccreting black holes with mass accretion rates 10 _1 — 10 Mq s _1 and 
viscosity parameters a vis — 0.01 — 0.1. Under such conditions the total neutrino luminosity emitted from 
the disk lies in the range L v = 10 51 - 10 54 erg s _1 . We focus our attention on the inner regions of the disk 
(within 10 Schwarzschild radii or so) where the midplane density exceeds about 10 9 gr cm~ 3 , and where the 
major fraction of the neutrino luminosity is generated. 

Now, for the range of temperatures considered the torus material should be a mixture of baryons and 
a light fluid (photons and electron-positron pairs in equilibrium). The light and baryonic fluids will be 
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tightly coupled in the sub-slow magnetosonic region, owing to the large Thomson depth. Deep beneath its 
surface, the torus is in a hydrostatic equilibrium where the vertical gravitational force exerted on it by the 
black hole is supported by the baryon pressure. In this region heating and cooling processes due to neutrino 
emission and absorption proceed at very high rates, and the temperature variation in the vertical direction 
is expected to be very small. As a consequence, the vertical density profile is roughly a Gaussian with a 
rather small scale height (h/R n £ 0.1; see below). As the density drops, the ratio of light fluid pressure 
to baryonic pressure increases roughly as T 3 /p 7 until at p 10 8 T 3 4eV gr cm~ 3 the light fluid pressure 
dominates over the baryonic pressure (see eq. [23] below). At this point hydrostatic equilibrium can no 
longer be sustained and the matter starts accelerating along magnetic field lines. This suggests that the 
mass flux carried by the wind should be controlled by the heating and cooling processes at the base of the 
wind. These processes are therefore incorporated into the flow equations derived below. To simplify the 
analysis we consider, in what follows, only time independent, axisymmetric models. As shown by Vlahakis 
& Koniglc (2003b), a relativistic MHD pulse can be adequately described by the same equations. 



2.1. Basic equations 

Since, as explained above, the light and baryonic fluids will be tightly coupled at the base of the 
outflow, they can be treated as a single fluid. We denote by p, p, e, h — (e + p)/pc 2 , the proper baryon rest 
mass density, total pressure, total energy density, and dimcnsionlcss specific enthalpy of this mixed fluid. 
The stress-energy tensor then takes the form: 

T af3 = hpc 2 u a U P +pg af3 + -^-(F aa F^ - -^-g af3 F 2 ), (1) 

where u a is the four-velocity measured in units of c, F^ v = d^A^ — d v A^ is the electromagnetic tensor, and 
g a P is the metric tensor. In what follows we take spacetime to have a Schwarzschild geometry, defined by 
the line element (is 2 = —a 2 dt 2 + aT 2 dr 2 + r 2 d0 2 + R 2 d(f) 2 , with a 2 = l — r s /r, and R~ rsinO. We denote by 
the source terms associated with energy and momentum transfer by some external agent. The dynamics 
of the MHD system is then governed by the following set of equations: energy and momentum equations 

-^={V~9T a(i ), a + T%T^ = qP, (2) 



continuity equation 



and Maxwell's equations, 



'-g 



(V^p U «), Q = 0, (3) 



F?~ = -^= g {V~9F f3a ), a = 4^/, (4) 

-Fa/3,7 + F&l,a + Ffa,0 = 0. (5) 

We assume the flow to be stationary and axisymmetric (d t =9^ = 0), and require infinite conductivity, 
u a F a p = 0. The above set of equations then admit two invariants: the mass-to-magnetic flux ratio 

n = <%- (6) 

Dp 

and the angular velocity, 

tl = 4*-- r ^- t F r ,=v*-%-%. (7) 
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Here u p is the poloidal 4-velocity of the fluid defined by u 2 = u r u r + u e ug, v p = Up/j is the corresponding 
3- velocity, where 7 = au 1 is the Lorentz factor, — /u t 1 B p is a rescaled poloidal magnetic field 
defined through B 2 — (B r B r + B e B e )/a 2 , and B$ — rF rd /a is the toroidal magnetic field. By contracting 
up with eq. (2), using the identity T^^upu^u 1 " = —upu u (uP)^, and the thermodynamic identity 
dh = dp/ pc 2 + (kBT/mis[C 2 )ds, where s is the dimensionless specific entropy and tun is the nucleon rest 
mass, we obtain the change in s along magnetic flux surfaces: 

(p/m N )k B Tu a s ta = -u a q a . (8) 

By contracting with eq. (2) and taking the t and <f> components, we obtain 

pc 2 u a £^ = -q t , (9) 
pc 2 u a C. a = (10) 

where 



^ = hu,-^F^ = hu,-—^, (12) 



and 

are the specific energy and angular momentum of the MHD system, respectively We henceforth assume 
that q<j, = 0, so that the specific angular momentum is conserved. By employing eqs (6) (7) (11) (12) we 
can express the toroidal magnetic field in terms of the flow parameters: 

47rr? a 2 C - R 2 n£ 
aR {a 2 — R 2 Q 2 — M 2 ) ' ( ' 

where M is the Alfven Mach number defined by M 2 = Anhrfc 2 / p = u 2 /u 2 A , with u\ = B 2 / (Airhpc 2 ). Using 
eqs (7), (11) and (12), and the normalization condition u a u a = —1 yields 

2 (£ nCfja 2 R 2 n 2 2M 2 ) [(C/R) 2 (£/a) 2 ]M* 

p h 2 (a 2 — R?Q, 2 — M 2 ) 2 ' [ ' 

We can differentiate eq. (14) along a given stream line to obtain 

(lnu p y = ^l, (15) 
where (') denotes derivative along the stream line \& =const, and 

D = -(a 2 R 2 n 2 M 2 ) 2 (u 2 p - u 2 SM ){u 2 p - u 2 FM )/u 2 A , (16) 
N = Ci(lnB p )' + C 2 (lna)' + &{\*R)\ (17) 
N q = U{]n£)' + &{\ns)'. (18) 

Here usm and ufm are the slow and fast magnetosonic wave speeds, respectively, and are given in eqs. 
(A18) and (A19), and c s is the sound 4-velocity defined by c 2 = a 2 /(I — a 2 ) with a 2 given in eq. (25) below. 
The coefficients Q are functions of the flow parameters, viz., Q = Q(Q, C,£, s, M), and are derived in the 
appendix. Equation (15) generalizes the result derived by Takahashi et al. (1990) to the non-adiabatic 
case. As seen, energy and momentum exchange with an external agent (i.e., N q ^ 0) formally modifies 
the conditions at the critical points. However, in practice we find this to be a small correction. As shown 
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below, the main effect of heating is to enhance the temperature and the specific entropy at the slow point. 
In cases where the stream function is known, the above set of equations, augmented by an equation of 
state h = h(p, s), can be solved to yield the structure of the flow. The poloidal magnetic field can then be 
expressed in terms of the stream function. Taking as the stream function ^(r, 9) we obtain, 

b p = • ( 19 ) 



2.2. Equation of state 

For the range of densities and temperatures considered here, we find that the relativistic particles in 
the vicinity of the slow magnetosonic point are non-degenerate. The total kinetic pressure of the mixed 
fluid is then given by p = pi + pb, where 

Pi = ^aT 4 = 1.2 x 10 26 T^ eV , dyn cm" 2 (20) 

is the light fluid pressure, and 

kT 

Pb = P W 26 p 9 T M eV dyn cm" 2 (21) 

m N 

is the pressure contributed by the baryons, with TueV being the temperature in MeV units, and pg the rest 
mass density in units of 10 9 g cm~ 3 . The dimensionless enthalpy per baryon of the mixed fluid is given by 



h = l + — + - — . 22) 

ptr 2 pc z 



It is convenient to define the thermodynamic quantity 

,,= ^=0.53^^. (23) 

Pb P9 

Clearly, the pressure is dominated by the light fluid when a > 4. In terms of a the change in entropy is 
related to the change in temperature and density through 

m N c 2 ( „ dp\ /3 o \ dT , ^dp 

It is readily seen from eqs (23) and (24) that a is roughly the entropy per baryon in regions where the 
light fluid pressure dominates; that is s ~ a in the limit a » 1. The relativistic sound speed can now be 
expressed as: 

2 _ / dln/A p b 5 + lQg + 2cr 2 
fls ~ [dlnp) s ~ pc 2 h 3(1 + 2a) ' ( > 

We also need the change in h at constant density. Using the above results we obtain 

dh\ kT 5 + 8o- 

~d~a) p ~ mArc 2 3(l + 2a)' ( ' 
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2.3. Rates for neutrino heating and cooling 

The neutrino flux emitted from the dense disk layers depends on the structure of and conditions in the 
disk, which are uncertain. Quite generally the neutrino luminosity comes predominantly from disk radii 
within a few r s (Popham et al. 1999). To simplify our calculations, we invoke a spherical neutrino source 
with a radius R v = 10 6 i?„6 cm and a total luminosity L v = 10 53 L„53 erg s _1 , which are treated as free 
parameters. Some fraction of the energy of the escaping neutrinos is deposited in the surface layers of the 
disk, via pair neutrino annihilation into electron - positron pairs and neutrino capture on neutrons and 
protons, thereby giving rise to a significant heating of the matter near the injection point of the wind. It 
is conceivable that turbulence or magnetic field dissipation may provide additional heating of the surface 
layers, but these processes are difficult to model and will be ignored. Detailed account of the various 
neutrino heating and cooling process is given in Qian & Woosley (1996). The heating rate due to neutrino 
absorption is approximately (see also Bethe and Wilson 1985) 

.3/2 

Q vn ~ 5 x 10 31 -if^p 9 / ra (r) erg cm" 3 s^ 1 , (27) 

and that due to neutrino-antineutrino annihilation into electron-positron pairs is (e.g., Goodman et al. 
1987; Qian & Woosley 1996) 

Q v - V ~ lO 31 ^/^) erg cm" 3 s" 1 . (28) 

Here f vn (r) = [1 - (1 - R 2 Jr 2 ) 1/2 }, and /„ v (r) = #„(r)[l - R 2 Jr 2 + 4(1 - R 2 Jr 2 )V 2 + 5]. Energy gain by 
neutrino capture on nucleons dominates at densities above 

p 9 = 0.2(L vB 3/ Rle) 1/2 (U/Un). (29) 

As a result of neutrino heating, the surface layers of the disk will quickly rise to temperatures in excess 
of several Mev, at which emission of secondary neutrinos by the inverse processes; electron and positron 
capture on nucleons 

e m ~ 10 27 p 9 Tl IeV erg s- 1 cm- 3 , (30) 
and electron-positron annihilation into neutrinos 

e vP = 5 x lQ 2i Tl eV ergs s" 1 cm" 3 , (31) 

becomes the dominant energy loss mechanism. Equating the last two rates one obtains the density below 
which the pair neutrino cooling rate exceeds the cooling rate due to electron and positron capture on 
nucleons, 

P9 = 5 x lO^T 3 ^. (32) 

The critical density is found to be always well above this value, implying that URCA cooling is the dominant 
energy loss mechanism in the regions of interest. The source term associated with energy transfer is given 

by 

(33) 

Well below the slow-magnetosonic point, where the flow velocity is small and adiabatic cooling is negligible, 
the temperature profile is determined the condition q° = 0. In this region the density exceeds the values 
given by eqs (29) and (32), and so the heating and cooling rates are dominated by capture (eqs. [27] and 
[30]). The temperature profile is then given to a good approximation by 



L 



1/4 



TW(r) = 6(-^) fV*. (34) 
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3. The Sub-Slow Magnetosonic Region - General Considerations 

3.1. The flow near the disk surface 

Deep beneath the slow magnetosonic point, the mater is close to hydrostatic equilibrium. Heating and 
cooling by the neutrinos proceeds at high rates and the temperature is maintained at the level given by eq. 
(34). If the pressure there is dominated by the baryons, then the density scale height at radius Rq is roughly 
h/R ~ {kTRo/GMmp) 1 / 2 ~ lO" 1 ' 5 ^/^) 1 / 2 ^/^, which is typically much smaller than unity near the 
black hole. Consequently, the base of the flow is located close to the disk midplane. Moreover, the ratio 
of the sound and Alfven speeds is c s /ua — 10^ 1 - 5 T 2 IeV / B p i 5 , so that under the conditions envisioned the 
Alfven Mach number is anticipated to be very small in the sub-slow region. To order 0(M 2 ) eqs. (All)- 
(A13) yield: 

Ci = -(a 2 -ii 2 n 2 ) 3 (l + u 2 )c 2 , (35) 
a 2 (a 2 ~R 2 n 2 )(£~n£) 2 
G = (l-a>)h* ' (36) 

Ca— C 2 ^. (37) 

Consider now a field line described by the equation R = R(z), and denote i?o = R(z = 0), R' — dR/dz\ z=0 , 
and N — (dz/dl)~ 1 N , where d/dl = (dz/dl)d/dz is the derivative along the field line, and ' denotes now 
differentiation with respect to z. Close to the disk midplane we can expand R(z) about z = 0. To first 
order in z eq. (17) yields 

N = Ci(InB p )' + ^ {(O 2 - n 2 )[R R' + (R' ) 2 z + R^z] + Q 2 [l - 3(R' ) 2 ]z} , (38) 

where fifc denotes the Keplerian angular velocity at R — i? , viz., Q 2 . — r s /2i?Q. At the base of the flow 
adiabatic losses are negligibly small and N q ~ 0. Now, in reality the poloidal velocity of the outflow should 
vanish at some height Zd(Ro) above the disk where the streamlines of inflowing matter joins the outflow. At 
this effective surface the angular velocity f2 is expected to be equal to the local Keplerian angular velocity. 
To leading order we then find fi 2 = £lf (1 — 3R' Zd/Ro), where the radius at which the field line meets the 
surface is related to Rq through Rd ~ R + R' Zd- Substituting the latter result into eq. (38) we obtain 

N = Ci(InB p )' + ^n 2 [z - 3(R' ) 2 (z - z d )\. (39) 

From eq. (16) it is evident that D < when u p < usm- Consequently, the condition N < must be fulfilled 
in the sub-slow region in order that the flow be accelerated. If the magnetic field lines are diverging, viz., 
(In £> p )' < 0, as one might expect, then the first term on the right hand side of eq. (39) is always positive, 
implying that the second term must be negative in the sub-slow region. As seen, there are two distinct 
cases. If the magnetic field lines are inclined at an angle larger than 7r/6 to the vertical (R' > l/\/3), then 
the second term on the right hand side of eq. (39) has a root at z\ = [3(i?Q) 2 /(3(i?p) 2 — l)]zd and it changes 
sign across. Along such field lines the outflow is centrifugally driven and can be initiated even in the cold 
fluid limit (at which (i = 0), as was first shown by Blandford & Payne (1982). The slow magnetosonic point 
is located close to the disk surface, at z sm < z\ (the exact location depends on the value of the sound speed 
at the disk surface). If the magnetic field lines are inclined at an angle smaller than ir/6 to the vertical, then 
the second term in eq. (39) is negative at z — and decreases linearly with z in the region where the above 
expansion holds. This corresponds to the regime of stable equilibrium, as defined in Blandford & Payne 
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(1982). Along such field lines the mass flux is thermally driven, similar to the case of a spherical wind. As 
shown below, the slow magnetosonic point in this case is located on larger scales, z sm <~ Ro, where the 
sound speed roughly equals the escape velocity. The above derivation is in accord with the results obtained 
by Ogilvie (1997) 2 . 



3.2. Conditions at the slow-magnetosonic point 

At the critical point D = and hence N + N q = 0. Neglecting terms of order M 2 and higher, one finds 



usm = c s and 



-/ 2 n l 



(1 + c 2 s )a 2 s (lnB P Y c + ^-(a c a' c - R c R' c n 2 ) - l c —^- = 0, (40) 

CK C ope Up 

where eqs (35)-(37) and (A14)-(A16) have been used. Henceforth, subscript c refers to quantities at the 
critical point. Equation (40) can be solved now for the sound speed a s . The result can be simplified 
further by noting that at the slow magnetosonic point the flow is sub-relativistic, so that c 2 ~ a 2 « 1, 
and 7 C ~ 1. Moreover, at the slow point the last term on the LHS of eq. (40) is approximately 
q t /3pc 3 u p ~ 10 _9 - 5 (L 1/5 3/i?^ 6 ) 3 / 2 (i? !y /r c ) 2 aJ 1 cm" 1 , where eq. (27) has been employed, and can be 
neglected to lowest order. With the above approximations the solution of eq. (40) reads: 

„2 _ „2 _ («c"c ~ RcR'SP) 

Now, in general we expect (lni3 p ) c = — A/R c , where A is an order unity number, with A ~ 2 for conical 
flux tubes, A = for cylindrical flux tubes, and A ~ 1 in the case of a self-similar magnetic field. In terms 
of ilk eq. (41) can be expressed as: 

_ r c f RcR'A^M - rtrp + n 2 z c Rl I ^ (42) 



where r 2 = R 2 ] + z 2 . From the last equation it is evident that a 2 « 1. Below we find that typically a s lies 
in the range 0.1 — 0.01. 



3.2.1. Centrifugally driven mass flow 

Since in this case the critical point occurs near the disk surface we can expand R c = R(z c ) about z = 
as before. Equation (42) then reduces to 

2 _ r s z c + 3(R' ) 2 (z d - z c ) 
asc ~Ro+r s AR ■ 

Adopting for illustration A = 2, r s /R n = 0.3, zj/Ro — 0.1, we estimate a 2 c < 10~ 2 at the critical point. 
Since the latter occurs close to the surface, adiabatic cooling is unimportant and the temperature profile is 
given to a good approximation by eq. (34) . For the range of neutrino luminosities considered here we expect 
T c £ 1 MeV. The mass loading of the outflow, p c a SCl would depend on the density at the flow injection 



2 Notc though that the term £i(ln_B p )' is neglected in eq. (5.35) of Ogilvie, which may be justified near the disk surface. 
However, the omission of this term precludes treatment of transonic outflows along low inclination field lines 
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point. However, the critical density and, hence, the mass flux cannot be arbitrarily small. To estimate the 
minimum mass flux, we employ eq. (25) to write 

, , 5 + 10<7 + 2er 2 

al = l^T MeV 3(i + 2a) , (44) 

and conclude that under the conditions envisioned a cannot be much larger than unity at the critical point 
in order for eq. (43) to be satisfied. Equation (23) then implies a critical density of p c £ 10 8 T ( ? Mey g cm~ 3 , 
and a mass flux M ~ p c a s irRQ > 10 30 g s _1 , for Rq = r s = 10 6 5 cm. Consequently, outflows along field 
lines having inclination angles larger than ir/6 are expected to be sub-relativistic, owing to the large mass 
flux driven from the surface. Numerical integration of the full equations confirms this result. 



3.2.2. Thermally driven mass flow 

In this case the critical point occurs higher above the disk. As will be shown in sec. 4, the specific 
entropy rises steeply during the initial acceleration phase, so that a >> 1 near the critical point. Using eqs. 
(25) and (41) we can relate the critical density to the critical temperature. In the limit a » 1 we obtain 

p c ~10 5 %2K g rcm -3. (45) 

ate 

As seen, the critical density is a sensitive function of the critical temperature and, therefore, adiabatic 
cooling can largely suppress the mass flux. An estimate of the mass-to-magnetic flux ratio can be obtained 
by noting that the critical point is located roughly at z c <~ Rq, so that B pc ~ B p0 . One then finds 

V= P -fr^ S^fsVfl- 1 gr cm- 2 s" 1 G _1 . (46) 

From the above result we further estimate that for the range of conditions considered here the relativistic 
energy per baryon of the matter is dominated by the nucleon rest mass, implying h c -f c a c ~ 1. To a good 
approximation the total energy per baryon (eq. [11]) is then £ = 1 + £g, with the magnetic energy per 
baryon given by 

E B ~ 5 x 10 2 Tj eV B 2 p015 a sc ^ . (47) 

At the highest neutrino luminosities adiabatic cooling can be neglected even along low inclination field 
lines, and the temperature at the slow point is given to a good approximation by eq. (34) with r = r c . The 
above quantities can then be estimated once the location of the slow point is known. At lower neutrino 
luminosities, however, adiabatic cooling will suppress the critical temperature, rendering such estimates 
highly uncertain. One must then integrate the flow equations in order to accurately compute the location 
of and temperature at the slow point. This is discussed in the next section. 

Equations (41) and (45) imply that, to order 0(M 2 ), the critical mass flux on a given field line, p c c Sl 
depends only on the gravitational potential, the angular velocity, the neutrino luminosity (through the 
critical temperature T C M e v), and geometry of streamlines, and is independent of the strength of poloidal 
and toroidal magnetic fields at the disk surface. Consequently, for a given choice of L u and magnetic field 
geometry, there is a critical value of BpoB^o above which Eg > 1 and the flow is essentially magnetically 
driven. Numerical integration of the flow equation confirms this conclusion, as shown below. 
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4. Numerical Model 

The model outlined above is characterized by three parameters: the black hole mass r s , the neutrino 
luminosity L u , and the neutrinospheric radius R v . Fixing these parameters, eqs. (8), (9) and (15) can be 
solved simultaneously once the magnetic field geometry is known. In general, however, the stream function 
is unknown a priori, and one must solve the highly non-linear trans-field equation simultaneously with eqs 
(8), (9), (22) and (24), and subject to appropriate boundary conditions to fully determine the flow structure. 
Unfortunately, the inclusion of gravity and external heating breaks scale-freeness, and precludes separation 
of variables as in self-similar treatments. This means that a self-consistent calculation of the field geometry 
is impractical. Since in this work we are merely interested in estimating the mass flux driven by neutrino 
heating it is sufficient to integrate the above flow equations only up to the slow magnetosonic point. As 
will be shown below, the latter is located rather close to the disk surface where the field geometry should 
anyhow reflect the boundary conditions. Our approach shall be to specify the stream function and use it 
to integrate the above set of flow equations along a given field line. We explore different field geometries in 
order to elucidate the dependence of baryon loading on the basic characteristics of the flux surfaces. The 
system under consideration has three integrals of motion: The angular velocity f2, the angular momentum 
C, and the mass-to-magnetic flux ratio, 77. In our approach, fl, C, and the values of £ and B p at the origin 
arc treated as additional free parameters, and i] as an eigenvalue of the system. 

In what follows we first describe the magnetic field geometries used in the our calculations. 



4.1. Topology of magnetic surfaces 

4-1.1. split monopole 

In our first example we consider a split monopole field, with its center shifted a distance a from the 
black hole along the hole rotation axis (see Fig 1). In cylindrical coordinates, the field line equation reads: 
R(z) = n(z + a) = kz + Rq, where n — tan# is the field line parameter and Rq = The stream function 
has the form \I/ = \& (1 — cos 9). Substituting the latter into eq. (19) yields the poloidal magnetic field: 



Bp 
~B~ 



a ^{(i+KV^(^+^r 3/2 } 



poj a 2 R 4 (l + K 2 )-r s /R 

where B p o = B p (z = 0). Note that this geometry implies a current sheet at the disk midplane. 
The proper length of a given stream line, measured in units of Rq, is given by 



(48) 



dl = ^(g rr dr 2 + gggdO 2 ) 1 / 2 = - 



a 



1 



r s R 2 



(R 2 + Z 2 ) 3 / 2 



1/2 



dz 

Rq ' 



and the derivative along stream lines by, 

a n Up d 
U da =R- dl 



r s Rl 



-1/2 



d_ 

dz ' 



(49) 



(50) 



{R 2 + z 2 f/ 2 _ 

From eqs. (9) and (50) we find that the change in specific energy per unit length along a given streamline is 



d£_ 

dl 



q t R a 
pc 3 u p ' 



(51) 
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Under the assumption that the neutrino trajectories are radial we obtain 



, , R dl\ 1 kR + z \ 

u vq =- q a{j-u p (a—) -j===\. (52) 



Substituting the last equation into eq. (8) yields 



ds q t R a \ ( dl\ 1 kR + z j 

7 — Up I a— I V , (53) 



dl kT(p/rriN)up j \ dz J ^/R 

which in the non-relativistic limit, 7 = a = 1, reduces to the result derived by Qian & Woosley (1996). 



4-1-2. Self-similar geometry 

In case of an outflow from a disk, we assume that the streamlines in the vicinity of some radius R = i? 
can be described by 

{R,z) = Ro\g(0,Z], (54) 

subject to the boundary condition = 0) = 1. The stream function is the solution of the equation 
rsm9/R (^) = g(r cos 6/ R ($>)). Using eq. (19) one finds 



B P V _ ocl f l + (.^) 5 



b p0 j k + (fto) 2 ]5 2 « 2 1 (9 - cs«) 2 i? (.9 2 + a 3/ 



(55) 



where subscript refers to the value of the corresponding quantity at the the disk midplane, viz., at £ = 0, 
and g^ — dg/d£. To obtain the derivative along streamlines we note that in terms of the proper length of a 
given streamline, R a dl = (g rr dr 2 + geedd 2 ) 1 ^ 2 , we have 

u a d a = u p d/dl = u p (dl/d£)- 1 d/d4, (56) 

with 



(w^li+ter-^-^^i^f. ,57) 



For the spherical neutrino source invoked above we obtain, 



^ = - q <* {l -u^y^0J>, (58) 



and 

d£ q t R a dl 
(ti; pc 3 u p 



(59) 



ds qtRo ) ( dl\ m±C{ , ar ., 

7 <*— - u p , , 5 ? • ( 60 ) 
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4.2. Boundary Conditions and Numerical Integration 

Equations (8), (9), (13), (15)-(18), (22)-(25), and (33) have been integrated numerically for different 
field geometries. In each run we fix the parameters B p0 , B<p , r s /R a , R v /R and L v53 /Rl 6 . In what 
follows, we find it more convenient to parametrize the neutrino flux in terms of the effective temperature 
T v ~ 6(L ly53 /i?2 6 ) 1 /4 MeV . The temperature at the origin (z = 0) is determined by the condition 
q t (z = 0) = 0, and is given to a very good approximation by To ~ T v fln (see eqs. [27] and [34]) . The 
integration starts at sufficiently dense disk layers where the pressure is marginally dominated by the baryon 
pressure. To be more concrete, for a given To the density at the origin po was chosen such that a £ 4. 
The entropy per baryon at the origin is taken to be s = 8.7 + \n{T^ v / p<j) + a (Popham et al. 1999). 
The parameter 77 is then adjusted iteratively by changing the boundary value of the poloidal velocity u p o, 
until a smooth transition across the slow magnctosonic point is obtained. To verify that the eigenvalue 
7] is insensitive to po 3 , each integration has been repeated several times for a given choice of our model 
parameters, each time with a different value of po- We find, indeed, that the dependence of r\ on p is very 
weak provided that the corresponding value of a is not too large. 



5. Results 

We now present solutions for the flow structure in the sub-slow region, for the two classes of field 
geometries described above. For the models calculated using the self-similar field geometry we invoked a 
parabolic field-line shape: 

ff(0 = l + ^ 2 . (61) 

where < S < 1 is an additional parameter, independent of Rq, that controls the divergence of magnetic 
field lines above the disk. Note that for this choice of g (In Bp)' = at z — 0. All the results presented 
below were computed using a black hole mass Mbh — 3Mq, and neutrinosphcric radius R v = 3r s = 10 6 5 
cm. For this choice of R v the effective temperature T v lies in the range 1 to 6 MeV for the range of 
neutrino luminosities discussed in Popham et al (1999). Sample results and the parameters employed are 
summarized in figs 2-8 and in table 1. 

Typical solutions are shown in fig 2, where the change along a field line of the flow quantities indicated 
is plotted against the normalized height above the disk midplane, z/Rq, for the split monopole geometry 
with k — 0.2, Ro = R u , _B p0 = 10 15 G, (—B < f,/B p )o — 0.1, and fi = fifc, and for different values of the 
neutrino luminosity given in terms of T v . For each of the cases shown the temperature at the origin is 
T = T„ since Rq — R v . As seen, the slow magnetosonic point is located rather close to the disk surface, 
at z ~ Ro, as one might expect. Well beneath the slow magnetosonic point the flow velocity is small, 
the density is high, and the specific entropy is dominated by the baryons, that is a 4. In this region 
the net energy deposition rate nearly vanishes, viz., q l ~ 0, and the temperature profile is given to a 
good approximation by equation (34). As the flow accelerates its temperature starts falling, and since the 
neutrino cooling rate is very sensitive to the temperature q t increases rapidly leading to a steep rise of the 
entropy per baryon to its terminal value, as seen in fig 2. The reason why the energy deposition per unit 
length peaks well below the slow point, as seen, is that u p there is small and, hence, the time over which 
a fluid clement is exposed to the external neutrino flux, dt ~ dl/u p , is large. Although formally the total 



3 The eigenvalue r) should depend, to some extent, on po through the condition q t (z = 0) = that fixes the the temperature 
at the disk layer where p = po • 
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energy per baryon is not conserved along streamlines, we find that it changes only slightly (to less than 
0.1 %) as the flow accelerates. This is because the net energy deposited per baryon comprises only a small 
fraction of the baryon rest mass. The main effect of neutrino heating is not to change the net outflow 
energy along streamlines, but rather to enhance the baryon load by increasing the light fluid pressure in 
layers of higher baryon density. Thus, in practice the total energy per baryon can be considered conserved 
in the relativistic and mildly relativistic cases (this may not be true if the wind is highly sub-relativistic), 
and so the value of £ presents, essentially, an upper limit for the asymptotic Lorentz factor of the wind. 
The values of £, r\ and the mass flux, which we somewhat arbitrarily define as M — pou p o2irRl (for a 
two-sided outflow), that corresponds to the cases depicted in fig. 2 are listed in table 1. The values of and 
£ computed numerically are in good agreement with analytic expressions given in eqs. (46) and (47). 

In fig. 3 we plotted solutions obtained for a centrifugally driven outflow along field lines in the 
unstable regime (k > l/y/S). The effective disk surface in this example is located at Zd = 0.1i?o- The 
temperature at the surface is T = 2 MeV, and the remaining parameters are the same as in fig. 2. As 
seen, the slow magnetosonic point is located very near the surface. The correspoding mass fluxes are 
M/(10 30 gr s _1 ) = 7.6, 76, 95, and 101 for k = 1, 2, 3, and 4, respectively. In addition to this example, 
we also explored other regimes of the parameter space. The main conclusion is that, quite generally, the 
mass loading of centrifugally driven winds from neutrino cooled disks is too large to allow acceleration to 
relativistic speeds. Since we are merely interested in relativistic winds here we focus, in what follows, on 
the regime k < 1/ y/3. 

Fig 4 exhibits the dependence of the total energy per baryon £ on the ratio {—B^/B p ) Q . It is clearly 
seen that the magnetic energy per baryon, £ B = £ — 1, depends linearly on (—B^/B^q, in agreement with 
the analytic expression given in eq. (47). The ratio of mass-to-magnetic flux rj is found to be practically 
independent of {—B^/B p )q. In the examples depicted in fig 4, rj changes by less than 3% over the range 
(—B t f,/Bp)o = 10~ 3 — 1. In fact we find that the mass flux pu p is highly insensitive also to the strength of 
the poloidal field B p o, provided the Alfven Mach number M is sufficiently small at the slow magnetosonic 
point (that is, c s « ua there), and that it depends predominantly on the neutrino flux, the geometry 
of magnetic field lines, and the angular velocity St. This confirms that for a given choice of the latter 
parameters, £b is proportional to the product {—B p B^)q. The results presented in fig 4 can then be readily 
rescaled for other choices of B p0 . 

The dependence of baryon loading on the inclination angle of the field line is examined in fig 5, where 
rj is plotted against the parameter n for a fixed i? (which corresponds to changing the distance a in fig 
1 keeping Rq fixed). The corresponding mass flux, M = pou p o2ttRq, is indicated on the right axis. As 
expected, the mass flux is larger on field lines with larger inclination angles, owing to the slingshot effect. 
The dependence of mass loading on the radius at which the field line intersects the disk is demonstrated 
in fig. 6, where rj and M are plotted against Ro for a fixed n. As seen, the trend is that the mass flux 
increases with increasing Ro. This mainly reflects the fact that the escape velocity is smaller on field lines 
that emanate from larger disk radii Rq. However, the dependence is weaker than one might anticipate, 
owing to the decrease in the heating rate with increasing Rq/Ru (see eq. [27] and below). We also made 
some runs with different values of r s jR v and found the mass flux to be sensitive also to this parameter at 
low T v . General relativistic effects appear to be important near the black hole, but we do not attempt here 
to quantify them. 

For completeness, we also examined the dependence of r\ on the angular velocity of magnetic flux tubes. 
Although the latter is naively expected to be nearly Keplerian, it is conceivable that various effects may 
give rise to non-Keplerian rotation of some flux tubes during certain periods. The general trend, as seen in 
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fig. 7, is for the baryon load to increase with increasing f2, by virtue of the larger contribution of centrifugal 
forces. However, a change in may also lead to a change in the Poynting flux, so the implications for the 
asymptotic Lorentz factor are not straightforward. 

Fig 8 presents results obtained using the self-similar configuration with <?(£) given by eq. (61), for 
S = 0.5 (left panels) and S = 0.2 (right panels). The behavior of the solutions is rather similar to the split 
monopole case. We generally find low mass fluxes for which £ >> 1 on field lines that emanate from within 
several r s and that diverge not too fast at low enough neutrino temperatures (T„ ^ 2 McV). We have also 
checked other forms of g(£). In all cases we find that at T v £ 2 MeV relativistic outflows can be launched 
from the innermost disk regions. 

We finally address the question whether the nuclear composition of the expelled matter changes during 
the acceleration of the flow. Significant evolution of the electron fraction Y e is expected if the neutronization 
timescale, t n ~ lOT^y s ' * s mucn shorter than the outflow time td = (Ro/c) J dl/u p . The evolution of 
Y e should take place at the base of the flow, where the temperature is highest and, hence, t n is shortest, 
and where u p is smallest and, hence, td is longest. For most cases studied above we find the outflow time 
to be comparable to t n , and so the n/p ratio is expected to evolve during the initial acceleration phase. 
Quantitative treatment of the nuclear processes in the wind is beyond the scope of this paper and is left for 
a future investigation. 

6. Conclusions 

We have explored stationary, axisymmetric GRMHD disk outflows in Schwarzschild geometry with 
neutrino-driven mass ejection, under conditions anticipated in hyperacrreting black holes. We have 
examined the dependence of mass loading of the outflow on the neutrino luminosity emitted by the disk 
and on other wind quantities, for a range of effective neutrino temperatures T v = 1 — 6 MeV, which for 
the simplified geometry of the neutrino emission zone invoked in our model corresponds to a neutrino 
luminosity range L v = 10 51 — 10 53 erg s _1 . We find that up to a few percent of the neutrino luminosity 
produced inside the disk is deposited at the upper disk layers and used up to eject baryons along magnetic 
field lines. The neutrino-driven mass flux depends predominantly on the effective neutrino temperature at 
the base of the flow, on magnetic field geometry in the sub-slow magnetosonic region, and on the angular 
velocity of magnetic flux surfaces, but is highly insensitive to the strength of poloidal and toroidal magnetic 
field components near the surface provided the Alfven Mach number is sufficiently small at the slow point; 
the dependence on neutrino luminosity and the shape of magnetic field lines near the disk surface appears 
to be particularly sensitive. The heating of the wind by the escaping neutrinos results in a steep rise of the 
entropy per baryon, s, during the initial acceleration phase, after which it saturates. The final value of s is 
larger for lower L v , but does not seem to reach extreme values. For the range of conditions explored above 
it is typically below lOOfes per baryon. 

Our principle conclusion is that ejection of relativistic outflows from the innermost disk radii, within 
several r s or so, is possible in principle for certain magnetic field configurations even in non-rotating black 
holes, provided the neutrino luminosity is sufficiently low, and the magnetic field is sufficiently strong . In 
the case of the split monopole field delineated in section 4.1, we obtain a mass flux of pu p ~ 10 14 gr cm~ 2 
s _1 on field lines in the vicinity of the innermost stable circular orbit that make an angle 9 <~ 12° with the 
rotation axis, for a neutrino luminosity L v ~ 10 52 erg s _1 ( corresponding to T„ ~ 2 MeV). The outflow 
along those field lines can, therefore, accelerate to Lorentz factors ;$ 100 if the magnetic field strength at 
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the disk surface is B p o ~ 10 15 G. The mass flux depends sensitively on both 6 and T v , as can be seen from 
figure 5 and table 1, and so a steep profile of the Lorentz factor across the polar jet is naively expected. 
We find a similar trend also for the other field geometries adopted. Thus, the picture often envisioned, of 
an ultra-relativistic core surrounded by a slower, baryon rich wind (e.g., Levinson & Eichler 1993, 2000; 
) seems a natural consequence of neutrino-assisted MHD disk outflows. At higher effective temperatures, 
T v > 2 MeV or so, the baryon load is typically well in excess of the value inferred in GRB fireballs, and the 
wind is expected to be slow. This conclusion may be altered if the central black hole is rapidly rotating. 

The hyperaccretion process is likely to be intermittent, leading to temporal changes in the neutrino 
luminosity. This should result in large variations of the Lorentz factors of consecutive fluid shells expelled 
from the disk in the polar region, owing to the sensitive dependence of mass loading on L v . In this 
situation we expect strong shocks to form in the outflow. If the polar disk outflow is associated with the 
GRB-producing jet, then the observed gamma-ray emission can be quite efficiently produced behind those 
shocks. 

It is worth noting that the conditions we find to be optimal for launching an ultrarelativistic jet 
in the polar region, arc also the conditions favorable for large neutron-to-proton ratio in the disk. For 
example, by employing the results presented in table 2 of Popham et al. (1999), we estimate T v < 2 MeV 
for accretion rates M acc ^ O.IMq s -1 into a 3 Mq Schwarzschild black hole. Under such conditions a 
neutron-to-proton ratio larger than 10 is expected inside the disk within a few r s for viscosity parameters 
a < 0.03 (Pruet et al. 2003; see also Beloborodoc 2003). If turbulent mixing can quickly lift up neutron 
rich matter to the surface layers, then in can be picked up by the outflow. Beloborodov (2003) has shown 
that deneutronization is not expected at temperatures below about 8 MeV if the mixing timescale is of the 
order of . However, because at low T v the flow time is comparable to the neutronization timescale, the 
electron fraction Y e may evolve as the flow accelerates and the final nuclear composition of the outflow may 
somewhat change. Further analysis is needed to determine the wind composition, but the trend seems to 
be that ultrarelativistic disk outflows can plausibly contain neutron rich matter. This is in accord with the 
conjecture often made, that fireballs pick up neutron rich material from the disk. 

I thank Arieh Konigle for inspiring conversations, the referee Nektarios Vlahakis for useful comments, 
and Jason Pruet for helpful correspondence. This work was supported by an ISF grant for the Israeli Center 
for High Energy Astrophysics 



A. Equation of Motion for General Relativistic, Neutrino- Assisted MHD Wind 

The projection of the momentum equation on the poloidal direction yields eq. (15). We can equivalently 
derive this equation by differentiating eq. (14) along a given streamline ^> = const (e.g., Camenzind 1986). 
To simplify the notation let us denote 

k Q = a 2 -R 2 n 2 , (Al) 
k 2 = (£-nC) 2 , (A2) 
C? £ 2 

(A3) 

Differentiating the above equations along a stream line, making use of the fact that rj, and C are conserved 
on magnetic surfaces, yields 

(fco)' = 2a 2 (lna)' - 2i? 2 ft 2 (lni?)', (A4) 
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(k 2 y = 2(£ - nc)£', 

2£ 2 C? £ 2 

(k 4 )' = - 2 %2 ( lni? )' - 2^ (lnf)'. 

a z R & a z 

Taking the enthalpy to be a function of the density p and the entropy per baryon s we further obtain, 

(ln/i)' = a 2 (lnp)' + |^(ln S )', 

(InM 2 )' = (In ft)' - (lnp)' = (a 2 - l)(lnp)' + |^(m S )', 

where a 2 = (dlnh/dlnp) s is the adiabatic sound speed. Upon differentiating eq. (14) along a stream line, 
using the above results together with eqs. (25) and (26), we arrive at 



(A5) 
(A6) 

(A7) 
(A8) 



(In «„)' = (N + N q )/D, 
where the terms in the nominator are given by 

N = Ci(InB p )' + C 2 (lna)' + &QnR)', 

with 



2\2 



Ci --(fco-M 2 ) 



(1 + Up)(k - M )c s - M 



,2 B l 



inhp 



1 f M e £ 2 
C2 ~ h 2 (l-a 2 ) 



- 3£ 2 - 



R 2 n 2 £ 2 _ 2C? 

~r? 



M 4 + a 2 fc 2 (3M : 



-*b)}. 



C 3 = 



/l 2 (l-a2) \ i?2 



£ 2 

-M 6 + ( 3£^ - 



a 2 C 2 _ 2R 2 il 2 £ 2 
~R 2 ^ 



M 4 -R 2 n z k 2 (3M 2 -ho) 



and 
with 



iV g = C4(hi£)' + C 5 (ms)', 







(fc - M 2 )[(fc - 2M 2 )(£ - fi£)£ + M 4 £ 2 /a 2 }, 
[k 4 M 6 - k 2 (k 2 - 3k a M 2 + 3M 4 )], 



^ 2 (l-a 2 ; 

(5 + 8a) sc 2 



(5 + 10cr + 2o- 2 ) h 2 



and the denominator by 

D = (a 2 - R 2 n 2 - M 2 ) 2 



(u 2 - c 2 s )(a 2 - R 2 il 2 - M 2 ) + M< 



B 2 



Anhp 



In terms of the slow and fast magnetosonic speeds, 

u% M = K- ^K 2 -c 2 u 2 A (a 2 -R 2 n 2 ), 



l FM 



K + JK 2 - c 2 u 2 A (a 2 - R 2 n 2 ), 



where 



K = - 
2 



(a 2 - R 2 n 2 )u 2 A + c 2 + 



B 2 



iirhp 



the denominator D can be rewritten as in cq. (16). 



(A9) 
(A10) 

(All) 

(A12) 
(A13) 

(A14) 

(A15) 
(A16) 

(A17) 

(A18) 
(A19) 

(A20) 
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Table 1: Numerical models 



Model a 


(k, S) h 


Ro/r s 


Ro/Rv 


L v 53/R v q 


T v 
MeV 


To 
MeV 


»? 

gr cm~ 2 s _1 G _1 


M 

10 30 gr s" 1 




£ 


SMI... 


0.2 


3 


1 


0.012 


2 


2 


0.18 


0.011 


44 


492 


SM2... 


0.2 


3 


1 


0.062 


3 


3 


7.4 


0.47 


32 


13 


SM3... 


0.2 


3 


1 


0.197 


4 


4 


116 


7.4 


29 


1.8 


SM4... 


0.1 


3 


1 


0.012 


2 


2 


0.028 


0.0018 


87 


3100 


SM5... 


0.1 


3 


1 


0.197 


4 


4 


22 


1.39 


55 


5.1 


SM6... 


0.3 


3 


1 


0.012 


2 


2 


1.6 


0.1 


19 


50 


SM7... 


0.3 


3 


1 


0.197 


4 


4 


1368 


86.6 


15 


1.03 


SM8... 


0.2 


9 


3 


0.012 


2 


1.25 


0.37 


0.024 


18 


161 


SM9... 


0.2 


9 


3 


0.197 


4 


2.5 


396 


25 


15 


1.1 


SSI... 


0.5 


3 


1 


0.012 


2 


2 


26.4 


1.67 


14 


4.3 


SS2... 


0.5 


3 


1 


0.062 


3 


3 


1008 


64 


13 


1.05 


SS3... 


0.2 


3 


1 


0.012 


2 


2 


0.29 


0.018 


29 


309 


SS4... 


0.2 


3 


1 


0.062 


3 


3 


313 


19.8 


13 


1.2 



a "SM" stands for split monopolc geometry, "SS" for self-similar geometry. 

6 The numerical values in this column correspond to values of k in the "SM" models, and to values of 5 in the "SS" models (see 
text for further details). 

c Soo refers to the asymptotic value of the dimcnsionlcss entropy per baryon s. 



-21 - 




Fig. 1. — Sketch of the split monopolc configuration employed for the numerical computations. The origin 
of the monopole is shifted a distance a from the black hole along the symmetry axis of the system. The 
magnetic field is supported by a current sheet in the equatorial plane. The field line parameter is defined as 
K = tan0. The radius Rq at which the field line meets the disk is indicated. 
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Fig. 2. — Profiles of various quantities in the sub-slow magnctosonic region, computed using the split 
monopole magnetic field of fig. 1 with n — 0.2, B p0 — 10 15 G, {—B^/ 'B p )q — 0.1, and Keplerian rotation, 
viz., Q = Ofe. Each panel corresponds to a run with a different neutrino luminosity (indicated in terms 
of the effective temperature T v \ see text for details). The quantities plotted in each panel are: the slow 
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Fig. 3. — Profiles of T/Tq and M$m for centrifugally driven outflows. The surface temperature in all cases 
shown is Tq = 2 MeV, and the remaining parameters are the same as in fig. 2. 
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Fig. 4. — Total energy per baryon, £, versus the ratio of the toroidal and poloidal components of the magnetic 
field at the disk surface, —(B < f > /B p ) - 
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Fig. 5. — Mass-to-magnetic flux ratio r\ (left axis) as a function of the split-monopole field line parameter k. 
The corresponding values of the mass flux M, as defined in the text, are indicated on the right axis. 
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B =10 G; (-B /B )=0.1; k = 0.2 

P o * p o 

10 r 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 



R7R 

v 

Fig. 6. — The dependence of r\ and M on the radius i?o at which the field line intersects the disk (see fig. 
1). The distance o in fig. 1 is adjusted such that n is kept fixed as Ro is varied 
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Fig. 7. — T] and M versus angular velocity of the field line, ft, given in units of the Keplerian angular velocity 
at the disk midplane Clk- 
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Fig. 8. — Same as fig. 2, but for a self-similar geometry with the parabolic field lines given in eq. (61). 



